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Electron dynamics in quantum dots manifests itself in spin-flip spectra through electric dipole 
spin resonance (EDSR). Near a neutrality point separating two different singlet charged states of a 
double quantum dot, charge dynamics inside a 2 x 2 singlet manifold can be described by a 1/2- 
pseudospin. In this region, charge dynamics is highly nonlinear and strongly influenced by flopping 
its soft pseudospin mode. As a result, the responses to external driving include first and second 
harmonics of the driving frequency and their Raman satellites shifted by the pseudospin frequency. 
In EDSR spectra of a spin-orbit couplet doublet dot, they manifest themselves as charge satellites 
of spin-flip transitions. The theory describes gross features of the anomalous half-frequency EDSR 
in spin blockade spectra [Laird et al., Semicond. Sci. Techol. 24, 064004 (2009)]. 
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Operation of electron spins in quantum dots is 
critical for spintronics 1 - and especially for quantum 
computing. 2-4 Electron spins can be operated by elec- 
tron spin resonance (ESR) driven by a time-dependent 
magnetic field B(t) 5 or by electric dipole spin reso- 
nance (EDSR) driven by a time-dependent electric field 
E(t)£r— As applied to nanostructures, EDSR is eas- 
ier to operate because it only requires voltages applied 
to gates rather than striplines and is usually more effi- 
cient than ESR. While the physical mechanisms of EDSR 
differ depending on the coupling of the electron spin 
s = icr to the electron momentum k or to its coordi- 
nate r through the spatially dependent Zeeman energy in 
inhomogeneous magnetic field B(r) (external, exchange, 
or hyperfine) jLiSril the frequencies ut of both ESR and 
EDSR are equal to the Zeeman energy ujb — \g\^BB/h 
due to the energy conservation. 

However, more recently Laird et al.— reported inten- 
sive EDSR in GaAs double quantum dots (DQD) at the 
half-frequency uj = \<jJb observed side by side with the 
EDSR at the basic frequency ui = lob^ it was detected 
by breaking Pauli blockade. They emphasized that a 
"half-frequency response is as far as we know unprece- 
dented in spin resonance" . Indeed, such a response re- 
quires strong nonlinearity that is unknown for spin reso- 
nance in its previous realizations. Also, such a behavior 
has no similarity in nonlinear spectroscopy where absorp- 
tion of the first and second harmonics is alternatively 
forbidden in centrosymmetric systems by parity conser- 
vation, a requirement reasonably well fulfilled in symmet- 
ric single quantum dots. 

Analysis of experimental datai 5 - implies a mechanism 
of this unique behavior. Most important, the half- 
frequency resonance was observed only in the spin block- 
ade region close to the degeneracy point of the singlet 
(1,1) and (0,2) charge states of DQD^dl This obser- 
vation suggests that the mechanism of nonlinearity was 
electrical and stemmed from the fact that the DQD was 
close to the resonance between two singlet configurations 



where the polarizability of the system was flight and 
strongly nonlinear. Also, the resonance was only detected 
in a device including a micromagnet. This observation 
indicates that the coupling of charge and spin dynamics 
underlying the EDSR was through r-dependent Zeeman 
energy. Indeed, the topography of stray field [Fig. 5(b) 
of Ref. 15] shows that B(r) changed considerably at the 
DQD-size scale. However, this specific mechanism of 
spin-orbit coupling is not essential for what follows, and 
in InAs where intrinsic spin-orbit coupling is two orders 
of magnitude stronger than in GaAs, similar effects are 
expected even in the absence of micromagnet. 

In what follows, I calculate single- and double- 
frequency electrical responses of a DQD in terms of elec- 
tric dipole moments r w (f) and r2 W (£) induced by a har- 
monic w-frequency perturbation. The result is specific 
for DQD's near the intersection of (1,1) and (0,2) en- 
ergy levels where the system is close to degeneracy; see 
Fig. 1. Away from this intersection, both dots react to 
the perturbation independently, and 

r(t) = -eE(t)/mcj 2 , (1) 

where ujo and m are the confinement frequency and elec- 
tron mass, respectively 8 ' 9 ' 19 Because of the large uq, the 
displacement of electron density is small compared with 
the dot size a w y^h/mojQ, nonlinear effects can be dis- 
regarded, and higher harmonics do not develop. 

The specific form of the charge stability diagram of a 
DQD in a wide parameter range depends on the shape of 
the dot, the potentials on the gates, and the electron- 
electron interaction inside the dot. However, as ap- 
plied to nonlinear responses that are the focus of the 
present Rapid Communication, only narrow regions near 
the charge-balance lines are of actual importance. In such 
narrow regions, all detailed parameters of the DQD can 
be absorbed into a few parameters of a model Hamilto- 
nian. To keep connection to the data of Ref. [l5|, I concen- 
trate on the vicinity of the charge balance line separating 
the stability regions of (0,2) and (1,1) singlet states that 
are displayed schematically in Fig. 1. In this figure, the 
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charge balance line is shown as a blue (heavy) line, and 
charge dynamics is considered in the direction perpen- 
dicular to this line. In the vicinity of the line, and not 
too close to its ends, the Hilbert space of the system is 
spanned by two charge states, (0,2) and (1,1). 

In the (0,2) and (1,1) basis, the Hamiltonian of a two- 
electron DQD is 



H 



-[e + h(t)]r z + vt x 



(2) 



where r is a vector of Pauli matrices acting in the singlet 
(0,2) - (1,1) subspace; the corresponding degree of free- 
dom will be termed a pseudospin in what follows. Here e 
is the static detuning equal to a half-difference of the en- 
ergies of (1,1) and (0,2) states, v is the tunneling matrix 
element, and 



h(t) = hcosut 



(3) 



is the time-dependent detuning. Because all dynamics 
proceeds inside the singlet subspace, the real spin is sup- 
pressed. 

Static part Hq of the Hamiltonian H can be diagonal- 
ized by an unitary transformation 



U = ar z + /3t x , a = V(A-e)/2A, = V(A + e)/2A 

(4) 

with the eigenvalues e± = ±A, A = Ve 2 + v 2 . There- 
fore, 2A is the excitation energy of the pseudospin degree 
of freedom. The sum J = A + e is usually considered as 
the exchange energy separating T Q triplet and the lower 
singlet S stated 

The transformed total Hamiltonian Hjj — UHU -1 is 



H v = At z + {eh(t)/A)T + (yh(t)/A)T x 



(5) 



where tq is a unit matrix in r space, and the dipole 
moment r(t) of DQD can be expressed in terms of the 
solution b(t) = (bi(t),b 2 (t)) T of a Schroedinger equation 



id t b{t) = Hub(t); 



(6) 



here and below U = 1. This dipole moment r(t) = r°(t) + 



' (t) consists of two terms, of which the first 

r°(t) - {a 2 \h(t)\ 2 + i3 2 \b 2 (t)\ 2 

+ a/9[6i(t)^(t) + b 2 (t)bl(t)]}Z 



(7) 



is large and scales with the interdot separation Z. The 
second term 

r s (t) = {2a(3(\b 1 \ 2 -\b 2 \ 2 ) 

- (a 2 - P 2 )[h(t)b* 2 (t) +b 2 (t)bl(t)}}Z s (8) 

scales with Zg = (Vtel^i + z 2 \ipu), where ^02 and ipu 
are wave functions of the singlet (0,2) and (1,1) states, re- 
spectively. This term is of the order of SZ, where S <C 1 
is the overlap integral of the single-electron wave func- 
tions ipL and ipR °f the left and right dots. The z direc- 
tion is chosen along the axis of DQD. 



Equation (J6)) cannot be solved analytically, and exper- 
imental data suggest a regime in which the frequency of 
the perturbation uj is large compared the energy scales 
of the static Hamiltonian, to ^> v,e. Indeed, this is 
a necessary condition for the energies of spin-blocking 
states (T+ , T_) be separated from unblocked singlet S 
and triplet To states by energies close to wb. Therefore, 
in what follows I perform an expansion in 1 [u and con- 
sider A/w, h/ui <C 1 as small parameters, while h and A 
may be of the same order of magnitude. 

To this end, the diagonal ±A terms in Eq. (JSJ can be 
eliminated by a substitution 



b 1 (t) = c 1 (t)e 



-iAt 



b 2 (t) = c 2 (t)e 



(9) 



Then nondiagonal terms of Hjj acquire factors 
exp(±2iAi), and all terms of the Hamiltonian become 
proportional to h. This transformation allows to take 
into account, in a consistent way, the slowly oscillating 
factors that become turn responsible for the Raman satel- 
lites. 

Next, both functions c m (t), m = {1,2}, are expanded 
in Fourier series 



n>0 



inujt 1 - 



f ), (10) 



where c mn depend on t but slowly compared with e ±ICJ *. 
Because the usual temperature of the thermostat ~ 100 
mK « 10/ieV is of the same order of magnitude as A, 
the state (cio,C2o) is a mixed one even when h — 0. It 
can be parameterized as c\q = cos</?, C20 = smLpe 1 ^, and 
this choice is general enough to provide results in terms 
of the density matrix approach. While the driving term 
h(t) results in a slow ^-dependence of (cio, C20), it can be 
disregarded in the calculations performed below. 

For calculating the response r(t) including the single- 
and double-frequency terms, dynamic equations for first 
two coefficients c mn (t), with n — 1 and n = 2, were solved 
by expanding in 1/ui. The solutions for c^ in include time- 
independent terms and terms oscillating as e ±2lA *. The 
appearance of the Raman satellites is seen most clearly 
from the explicit form of the Fourier components 6^ 
cf x e- iAt and bf t ' " 



c 21 e 



iAt 



b± 



T—r—e ±——e e^smifi, 
2Aw 2Aw ^' 



v h 
2Acj 



iAt 



±-L^e iAt e^sin^ 
2Aw 



(11) 



that include both e ±lA * harmonics. Next, the Fourier 
components, bf 2 and b\\-, can be expressed in a simi- 
lar way in terms of 6^ and b\\ with corresponding time 
dependent coefficients. Expressions for the uj- and Iuj- 
responses of Eqs. (J7J and (|SJ| are more cumbersome be- 
cause they include biquadratic forms of the amplitudes 
bi{t) and b 2 (t). 

Before discussing specific results for these responses, it 
is instructive to take into account that r (t) satisfies, due 
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to the equations of motion of Eq. (J6)), an exact identity 

id t r°(t) = v(b 1 b* 2 -b* 1 b 2 )Z. (12) 

The derivative dtr°(t) is proportional to u> both for the 
first and second harmonic responses. Therefore, the r°(t) 
term in r(i), while scaled with a large length Z, includes 
an additional w factor in the denominator and therefore 
decreases with ui faster than r s (t). 

The w-responses calculated using Eqs. ([7]) and ((U, in 
the leading order in 1/uj, are 

r uj(t) — H — 'oi v cos 2y + esin2(^cos(0 + 2 At)]Z cos ut, 
Alu z 

(13) 

r%(t) = -— sm2cpsm((f) + 2A£)Zssxnut. (14) 

CO 

Therefore, the linear response includes oscillations at 
the driving frequency lj and two satellites at frequen- 
cies ui± 2A. The satellites originate from slowly oscillat- 
ing factors e lAt in the Fourier amplitudes of the same 
type as in Eq. flTTl) . To evaluate the relative magni- 
tude of both contribution, one can estimate v ~ Ssq, 
where eo ~ 1 meV is the intradot electron energy. Then 
r tl r t ~ <^b/ £ o- With uj b ~ 2.5 /ieV for GaAs at 
B ~ 100 mT, r° appears as a dominating contribution 
despite its faster decrease with u>. 

Double frequency responses calculated in a similar way 
are 

h 2 v 

r L(*) = 3 sm2ipsin((f> + 2At)Zsin2ujt, (15) 



r 2i» (*) — X 2 [ u cos 2^ + e sin 2<p cos(0 + 2 At)] Zs cos 2cj£ . 

(16) 

For the reasons explained above, the r 2ul term is expected 
to dominate, and it does not include the unshifted (cen- 
tral) second-harmonic contribution but only a doublet of 
Raman terms with the frequencies 2(w± A). Meanwhile, 
r 2bJ includes a central 2w line and its 2{ui ± A) satellites. 

Notably, 2 A energy shifts of the oscillation frequencies 
are universal for all satellite bands, both in the linear 
and quadratic responses. Because 2A is the energy sepa- 
ration between e± sublevels of the static Hamiltonian Hq, 
the origin of these shifts can be easily understood. Phys- 
ically, they come from the flopping between two pseu- 
dospin levels driven by the high-frequency perturbation 
h(t). Technically, they originate in the biquadratic ex- 
pressions of Eq. ([7]) and © from multiplication of e ±iAt 
factors in the amplitudes b(t). 

A comparison of Eqs. (fTB"]) and (fl"4")l and Eqs. (fT5")) and 
([T^l shows that the ratio r 2w /r u is of the order of h/oj, 
i.e., it is small in the parameter h/ui used above for de- 
riving all equations. However, because h appears only in 
prefactors rather than in the functional dependencies, it 
is probable that all qualitative conclusions are valid until 



h ;$ co. While the magnitude of the third harmonic has 
not been calculated explicitly, its magnitude (and there- 
fore the intensity of the uj = \ojb resonance) is expected 
to include an additional power of h/to. 

Coexistence in a geometrically symmetrical DQD of 
the first and second harmonics stems from the fact that 
electronic density near the charge balance point of (0,2) 
and (1,1) configurations is asymmetrical. While the elas- 
tic term in the nonlinear response of Eq. (|16[) may be, 
at least partially, attributed to the static asymmetry of 
electron density distribution, the Raman shifted terms 
originate completely from charge dynamics in the pseu- 
dospin space. Dynamical asymmetry manifests itself ex- 
plicitly in strong parity violation of 2A satellites that are 
ubiquitous and appear both in the w and to 2 responses. 

All dipole moments in Eqs. (fi~3|) and (fl6|) vanish after 
averaging over tp. This means that their amplitudes and 
phases strongly depend on the quantum state (cio,c 2 o)- 
Whenever the dynamics driven by the field h(t) is faster 
than the change of this state controlled by its coupling 
to the bath, both linear and nonlinear polarizabilities 
remain unchanged. On a longer time scale the polariz- 
abilities fluctuate and may be estimated by their mean 
square values. 

Now we are in a position to make connection to the 
data of Ref. HH. According to the above assumptions, 
oscillations of electronic density parameterized above in 
terms of the dipole moments r(t) drive EDSR through 
the B(r) mechanism and in this way affect the efficiency 
of the Pauli blockade. In Fig. 6 of Ref. [HI both the lin- 
ear and quadratic responses are split, at least into two 
components (visibility of the components is reduced by 
fluctuations of the nuclear bath). At B 100 mT, where 
this splitting is best seen, it is ~ 10% of the Zeeman en- 
ergy. Because the splitting does not increase with the 
external magnetic field B, it can not be related to the 
difference in the ^-factors of two dots. Therefore, it is 
natural to attribute it to the satellite structure unveiled 
above, which physically corresponds to two charge (pseu- 
dospin) Raman-type satellites accompanying each spin- 
flip line. Estimating the splitting as 4A in the spirit of 
Eqs. (fT5|) - (TTB jl . we arrive at A ~ 0.1/ieV which corre- 
sponds to tiny tunnel matrix element v and detuning e. 
The visibility of both resonances, especially of the half- 
frequency lu = \ojb resonance, decreases with j5, in a 
qualitative agreement with the decrease of the dipole- 
moment matrix elements of Eqs. (I13[) - (I16I) . 

While the gross features of the exotic EDSR observed 
in Ref. [TH are in a reasonable agreement with the above 
theory, the real situation is more complicated. Indeed, 
while the lu — ljb resonance manifests itself as a peak 
in the current, the u = |ub resonance is seen as a 
dip in the background current leaking through the Pauli 
blocked DQD. Most probably, this indicates that the 
half-frequency feature appears as a Fano resonance,^ 
and therefore different mechanisms of breaking the spin 
blockade (that is weak in the vicinity of the charge de- 
generacy point) should be taken into account. Moreover, 
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signs of both resonances change depending on the ex- 
act gate configuration,— which is not surprising having 
in mind the tiny magnitude of A estimated above. De- 
tailed patterns of the Pauli blocked currents near the 
charge balance point are outside the scope of the current 
paper. 

The mechanism of EDSR described above allows for (i) 
a new spectral region of the operation of singlet-triplet 
DQD qubits and (ii) a strong and controlled dynamic 
coupling between them. 

In conclusion, near the charge neutrality point the 
polarizability of a double quantum dot is large, highly 
nonlinear, and strongly influenced by dynamics of its 
soft pseudospin mode. A high-frequency detuning sig- 
nal h(t) = h cos ujt drives dynamics inside the 2x2 sin- 
glet subspace producing flops of the pseudospin mode be- 
tween its two eigenstates separated by energy 2 A -c oj- 
This charge dynamics results in oscillations of electron 
density at the frequencies w and 2oj and at the Raman 
shifted frequencies w±2A and 2w±2A. Pseudospin flops 



also act as an efficient mechanism of parity violation. 
Through one of the mechanisms of spin-orbit coupling, 
the density oscillations drive EDSR, and shifted frequen- 
cies manifest themselves as Raman-type satellites reveal- 
ing the pseudospin charge dynamics in spin-flip spectra. 
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Caption to Fig. 1 

Fig. 1. (Color online) Charge-balance diagram in the 
vicinity of the (1,1) - (0,2) charge-balance line shown in 
blue (heavy solid). Both time-independent e and time- 
dependent h(t) detunings are applied along the horizontal 
dashed line marked as e. Along this line, charge dynamics 
within the 2x2 (l,l)-(0,2) singlet subspace is described 
in the framework of a pseudospin- 1/2 formalism. 
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